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TECHNICAL NOTE D-685 


AN EXPLICIT LINEAR FILTERING SOLUTION FOR 
THE OPTIMIZATION OF GUIDANCE SYSTEMS 
WITH STATISTICAL INPUTS 
By Elwood C. Stewart 

SUMMARY 


The determination of optimum filtering characteristics for guidance 
system design is generally a tedious process which cannot usually he 
carried out in general terms. In this report a simple explicit solution 
is given which is applicable to many different types of problems. It is 
shown to be applicable to problems which involve optimization of constant- 
coefficient guidance systems and time-varying homing type systems for 
several stationary and nonstationary inputs. The solution is also appli- 
cable to off-design performance, that is, the evaluation of system per- 
formance for inputs for which the system was not specifically optimized. 

The solution is given in generalized form in terms of the minimum theo- 
retical error, the optimum transfer functions, and the optimum transient 
response. The effects of input signal, contaminating noise, and limita- 
tions on the response are included. From the results given, it is possible 
in an interception problem, for example, to rapidly assess the effects on 
minimum theoretical error of such factors as target noise and missile 
acceleration. It is also possible to answer important questions regarding 
the effect of type of target maneuver on optimum performance. 


INTRODUCTION 


There are a number of statistical theories which have been derived 
in recent years for the purpose of determining optimum system design and 
optimum theoretical performance. Three problems are invariably encountered 
in the application of these theories. The first problem is that the solu- 
tions generally involve long and tedious computations. Moreover, explicit 
solutions can rarely be obtained because of the complexity of the equations 
a.nri the many factors involved. For this reason solutions are usually 
carried out numerically for specific cases. Such a procedure contributes 
little to a basic understanding of the problem and makes it difficult to 
draw general conclusions as to the relation of the guidance and control 
task to both the best theoretical performance which can be achieved and 
to the optimum system design. Such considerations are important, for 
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example, in evaluating the effect on theoretical minimum error of such 
factors as limited vehicle maneuverability, the amount of noise, etc. 

A second problem is concerned with off -design performance, that is, 
the evaluation of system performance for inputs for which the system was 
not specifically designed. In general the actual input to a system will 
be different than the design input. For example, the actual noise level 
might be different than design noise level, or the actual and design sig- 
nal inputs might even be different processes. For this and other cases 
it is important to evaluate the deterioration in error for these off-design 
conditions . 


A third problem in system optimization concerns the choice of signal 
component of the input, which in the interception problem is the target A 

motion. Since the target may move in many different ways, it is important 2 

to consider two aspects: (l) the effects of type of signal characteristic 7 

for which the system is optimized, and (2) the effect on performance of 1+ 

subjecting the system to signal inputs for which the system was not 
specifically designed. 


This report will be concerned with the above three problems. The 
first two sections are concerned principally with the first problem. In 
the first section will be derived explicitly a simple but approximate 
solution to the filtering problem. Because of the length of this section 
a resume is given at the end. In the second section will be shown the 
applicability of this solution to many different problems of interest; 
such problems involve the optimization of constant-coefficient guidance 
systems and time-varying homing systems, for stationary and nonstationary 
signal inputs. The third section is concerned with the off-design problem 
(the second problem) . The last section is an example concerned with the 
third problem, the effect of type of signal input on the performance of 
optimum and nonoptimum systems. 

LIST OF IMPORTANT SYMBOLS 

a,p acceleration of target, ft/sec 2 

A 2 mean- square acceleration, ft/ sec 2 

c steady- state output defined in equation ( 4 ) 

o 2 

E mean- square error , ft 

optimum compensating network transfer function 
Hf fixed network transfer function 

quantity related to c^ (see eqs. (b62) ) 

N noise magnitude or zero frequency spectral density, ft 2 / radian/ sec 


< OJ C—-4- 
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p 

weighting function 

Q 

quantity related to c G (see eqs. (B 2 l) and (B24) ) 

\ r 

restricted quantity 

R 2 

mean- square restricted quantity 

s 

variable in the Laplace transform 

T 

A 

time at collision, sec 

/i 

2 ^o 

unit impulse 

7 

It- x 

dimensionless frequency, u/p (see eq. ( 28 )) 

y 

output in adjoint diagram sketch (c) 

• Y 

x O 

optimum closed- loop transfer function 

a 

constant multiplier in cp (see eq. ( 23 )) 

3 

input parameter, l/sec (see eq. ( 28 )) 

7 

vehicle parameter, sec (see eq. ( 37 )) 

€ 

error, ft 

n 

dimensionless parameter in optimum transfer function, 7 P 

*k 

the kth determinant 

V 

input dimensionless parameter, |/p 

i 

input parameter in frequency characteristic of i|r s , l/sec (see 
eq. ( 21 )) 

p 

Lagrangian multiplier, sec 4 

a 

constant multiplier in \|r s , ft 2 / sec 3 (see eq. ( 2 l)) 

9 

frequency factor of input 

♦ 

frequency factor of input analogous to spectral density 

• 

U) 

angular frequency, radians/ sec 

( ) 

complex conjugate of ( ) 

( ) 

ensemble average of ( ) 


Subscripts 


s 

signal 

r 

restricted quantity 

a 

acceleration 

n 

noise 

G 

error 


EXPLICIT FILTERING SOLUTION 
System Description 


The type of guidance system with which we will "be concerned has been 
described previously in references 1 and 2, but for completeness a brief 
summary will be desirable here. A great many guidance and interception 
problems can be represented by the block diagram shown in sketch (a) . 

This representation can be shown to apply, for example, to certain space 



navigation and missile problems which are time-invariant, or even to time- 
varying navigation problems, as will be shown later. It does not apply 
however, to guidance systems which operate part of the time as an open- 
loop system, such as certain interception weapons systems, fire control 
systems, etc. The two inputs to the guidance system are the signal x s , 
which contains the true information about the motion of the target to be 
intercepted, and the noise x n , which enters unavoidably with the desired 
signal. The outputs of interest are two in number. The first is the 
error e which is a measure of how far the system output deviates from 
the desired signal part of the input. The second is the saturating quan- 
tity r. Although there may be many quantities subject to saturation in 
a guidance system, it has been shown (see ref. l) that only the most 
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critical quantity need "be considered. In the interception problem this 
quantity is the thrust producing quantity, such as the deflection of a 
gimbaled rocket engine or of a control surface of a missile. In order 
that this quantity will appear explicitly, the guidance system is split 
into a compensating filter and an output element. The output element, 
which represents the dynamics of the vehicle, is assumed to he the fixed 
element because it is relatively unalterable compared to the remainder of 
the system. The other element is the compensating network which the 
designer is free to choose. 

The above representation for the guidance system does not represent 
the actual physical form of any particular system; that is, any linear 
missile guidance system of the type discussed above regardless of its form 
A or the number of feedback loops can be put into the equivalent form. shown. 

2 Conversely, for any given characteristics of the compensating and fixed 

7 filters, there are any number of corresponding physical systems which can 

4 be constructed. 


Performance Equations 


The purpose of this section will be to present the exact performance 
equations, that is, the expressions for the error and restricted quantities 
which will apply for any given system. These equations will serve as a 
starting point for optimization problems as well as off -design problems. 

The diagram shown in sketch (a) can always be redrawn into the form 
shown in the following sketch. Here the two inputs are taken to be 



Sketch (b) 


impulsive in character. Such a representation is a natural one both 
because nearly all inputs can be represented by a series of impulses and 
because the adjoint theory to be used later requires impulse type inputs. 
Each of these two inputs is then modified by a frequency factor cp(.sj 
before entering the system transfer function w(s) . The actual signal 
input, x s (t) , to W s (s) maybe an analytic signal, such as a step, etc., 
or some kind of random process. However, the noise input, x n (.t; , to w n (.s; 
will be restricted to be a stationary random process, which occurs so 
often in practice. The output, 0 o (t) , may represent either. of the two 
outputs of interest, the error or the restricted quantity , if W s and W n 

are properly defined. The mean- square ensemble average, 0 O 2 (T) , at a 
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particular time T is given "by the following expression, assuming no 
correlation between the two inputs: 


0 O 2 (T) = 0 S 2 (T) + © n 2 (T) 


( 1 ) 


Each of these components can be evaluated by use of the corresponding 
adjoint diagram (see refs. 3 and 4 for a description of the adjoint theory) . 
The adjoint system corresponding to sketch (b) is shown below where it can 
be seen that now there are two outputs, y n and y s , and only one input, Uq, 


y n (T) 


y s ( T > 



Sketch (c) 


JL 


which is always an impulse . For the moment let us confine our interest 
to only one of these components, the kth. The mean-squared ensemble 
value of the output 9k(t) at a particular time T can be expressed in 
terras of the adjoint response y^r) by the following general expression 


0 k 2 (T) 


yk 2 ( T ) p k(r) dT 


( 2 ) 


where P k (T) is a weighting function dependent on the nature of the input 
process in the real time domain. Such an expression is valid for linear 
systems which are subjected to either analytical inputs, such as a step, 
impulse, etc., or to random processes. For example, for a stationary 
random process, P(r) = while for a step input which starts at ti, 

P(t) = Uo[ T -(T-ti) ] , the -unit impulse. For the situation illustrated in 
sketch (c) , we may now use the relation (2) to rewrite (l) as 


0 o 2 (t) = 


vT pT 

ys 2 ( T ) p s( T ) dT + / y n 2 ( T ) p n( T ) dT 


(3) 


It often happens that because of the nature of the input frequency func- 
tion cp s (s) and the system W s (s) , the output y s ( T ) has a steady-state 
value c . Thus it will be found convenient to let 


y t ( T ) = y s ( T ) - c 


W 


1 The 2n is due to definitions of spectral density (see ref. 4, p. 9) • 
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A 

2 

7 

4 


Since in the cases in which we will be interested P^(t) is independent 
of T , equation (3) becomes 


0 O 2 (T) = 


yt 


'(t) dT 


2P s c 


yt( T ) dT 


+ P c c 2 T 


'T 


+ p n / y n 2 ( T ) dT (5) 


By the definition in equation (4) , we know that yj. -*• 0 as T -*• oo. Further- 
more^ in the usual case , y( t) becomes so small for reasonable values of T 

pOO p T 

that / y t 2 (r) dr « / y^ 2 (T)dT. Physically, this inequality corresponds 

J T 'Ad 

to the assumption that the system response time is less than the interval 
of interest T. Thus by increasing the upper limits to oo and by trans- 
forming to the frequency domain 2 , we can re-express equation ( 5 ) as 


0 o 2 ( t ) 



Y t (“) 


2 dw + 2P s c lim Y-^(s) 
s-»o 


+ P S c 2 T 


+ 



cp n ( w ) w n( w ) 



(6) 


An expression for Y-^(s) can be obtained from equation (4) and sketch (c) : 


where 


Y t (s) = cp s (s)W s (s) - § 

c = lim scp s (s)w s (s) 
s-^o 


Thus, 


(7) 

(8) 



c 

ioo 


2 

dw + 2P s c lim 
s->o 


Ts( s)W s ( s ) 


c 

s 


+ P S C‘ 


+ £n 

2jt 



cp n (w)W n (w) 


2 

dio 


(9) 


Equation ( 9 ) gives the mean-square ensemble or time average of any quantity 
at time T in terms of the system transfer function and the input weight- 
ing and frequency functions. 


"Note that 2 tc / |y(T) | 2 dr = / | Y(w) | 2 doj. 
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Now the expressions for the error e and restricted quantity r 
can he written. We introduce E 2 and R 2 to represent either mean-square 
time or ensemble averages of the error and restricted quantity, respec- 
tively. If sketch (a) is interpreted in terms of sketch (b) , we can 
deduce the following. 


For error € : 


For restricted quantity r: 


W s (s) = l-Y(s) 
W n (s) = Y(s) 

W s (s) = Y(s)/H f (s) 
W n (s) = Y(s)/H f (s) 


( 10 ) 


Therefore the error and restricted quantities become, from equation ( 9 ) > 


* - K 


cp s (w)[l-Y(w)] - ^ 


dw + 2P s c e lim ^cp s (s) [l-Y(s) ] - 
s-*o 


p p 00 

+ P s c e 2 T + ^ / |cp n (w)Y(w) | dw 

-.00 


( 11 ) 




°0 

cp s (w)Y(w) 

: L 


Hf(w) 

Pn r 

X) 

cpn( w ) Y ( w ) 

2 *L 


H f (w) 


10 ) 


dw + 2P s c r lim 


s-*o 


(p s (s)Y(s) 
Hf(s) 


^ !• + P S c r 2 T 


dw 


( 12 ) 


A 

2 

7 

k 


These two expressions will he useful throughout the remainder of the 
report . 


Optimum System Transfer Function 


In this section will he given a simplified explicit solution for the 
optimum system based on certain simplifying approximations. 

Exa ct equations .- The optimization problem is one of determining 
the optimum over-all system transfer functions which will minimize 

E 2 + pR 2 


(13) 


•£-->3 ro > 
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where p is a Lagrangian multiplier having the units of (e/r) 2 , and E 
and R are given by equations (ll) and (l2) . Such minimization could 
probably he carried out by brute force. However, anticipating the nature 
of the final solution can simplify the problem. After a little reflection 
one might expect that both the steady-state terms, c e and c r , would be 
zero. We will assume that this is the case and after we have the answer, 
we can check back to verify the assumption. Equations (ll) and (12) now 
become 


E 2 = 


1-Y(u>) I 2 


P S I tPsCkQ 1 

2rt 


- 00 


dm+ / | Y(co) | 2 P - n - M u) r 


dco 


(lb) 



P si 9s( tAj ) I 

2j t 


du + 


pOO 

Y(w) 

2 Pnl9n(w) I 2 

L 

%(<*>) 

2jt 


dco 


( 15 ) 


The first term of each of these expressions is due entirely to signal, 
the second to noise. Although these equations are written in terms of 
the fixed filter Hf(w) and the over-all system transfer function Y(go) 
(instead of the compensating network as described in sketch (a)), they can 
be expressed in terms of the fixed filter and the compensating filter by 
the relation 


Y(cj) = H c (o))H f (w) 


(16) 


An expression for the optimum compensating network which minimizes 
(13) can be readily found. If we make the following definitions 


t s (w) 

t n (w) 


P s l ^Ps( w ) I 

2* 

P n l c P n ( aj ) I 2 

2it 


( 17 ) 

(18) 


then it is clear that the \|r* s are analogous to spectral densities and 
that the solution for the optimum compensating network function, 1 ^ 0 ( 10 ) , 
is given by (see ref. 5 or l) 


Hco(iw) 


1 r -iwt r H f (ia)t s (a)e lat 

2jtA + (w) J o S A" (a) 


da dt 


( 19 ) 


where 


A (w) = [H f (iw)H f (iw) + p][\jr s (co) + t n (w) ] 
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As discussed in reference 1, the optimum compensating network Hco 
only he determined "by the simultaneous solution of equations (15) and (19) 
so as to result in the desired restriction. Evaluation of equation (l4) 
then follows. Of the various factors involved in these equations the 
quantities which must normally "be known are the input quantities \|r s (w) 
and 4r n (w) , and the fixed network Hf(to) . Unfortunately even for very 
simple forms for these functions, the complexity of operations involved 
in the solution of these equations does not permit the general solution 
to "be obtained explicitly. We will "be interested in certain forms for 
these functions which will now "be discussed. 

Simplified forms for inputs and fixed network .- A very common and 
important form for the inputs which occurs in many physical problems and 
particularly in the interception of targets is the following 


*sH = 2 ,2 

(20) 

W + I 


+s( W ) = 4, 2 ,2v 

CO (u) +£ ) 

(21) 

tnM = N 

(22) 


where tyg(u)) corresponds to the second derivative of the signal. The form 
in (20) or (21) is valid for many stationary and nonstationary processes 
as will be seen. Furthermore this form is general enough to approximate 
a variety of experimentally determined input data. The other function 
^n(co) , the noise frequency factor, is approximated by a constant. This 
is a good approximation, since in most physical situations the bandwidth 
of the actual noise is much broader than that of the optimum system trans- 
fer function Y 0 (w) . Obviously, from equation (17) the actual filters 
which appear in sketch (c) would be 

cp s (s) =i/gH/h = a.f(s) (23) 

s 2 (s H) 


cpn(s) = n/1 


'Phe functions cp and are merely two different ways of expressing the 
signal input function. Thus the gain a of the actual filter in 
sketch (c) is related to the magnitude a of the input frequency function 


a 


2rt 


(25) 
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A 

2 

7 

h 


The other function, Hp, is the fixed element representing the dynamics 
of the vehicle. This function is generally very complicated; when the 
output is taken to he a displacement quantity as it is here, the fixed 
network will he of the form 


%(s) 


= kj; 


a n s n +a n _ 1 s n_1 + . 
s 2 ( b m s m +b m .. 1 s 131 " 1 + 


• +1 

. • +l) 


( 26 ) 


Even for cases in which m and n are only 1 or 2 , the solution of equa- 
tion (19) cannot he achieved in general terms. However, there are several 
sound reasons to believe that such a complicated Hp is hoth unnecessary 
and undesirable, and that it can he approximated by 



( 27 ) 


The reasoning is as follows. First, it can he shown that the fixed net- 
work given by equation (27) is optimum. In this sense it should he noted 
that we are considering optimization within a class of opt imu m systems. 

That is, for each fixed network Hp there is an optimum over-all guidance 
system and corresponding minimum error. However, the relationship between 
the fixed network and the minimum error is so complicated there is no a 
priori way to tell which fixed network will he best. Nevertheless, the 
results of reference 2 show that the Hp given by equation (27) results 
in an over-all system which is the best of all these systems. Thus (27) 
is a desirable form for the fixed network and should be striven for. 

Second, even if the dynamics are not the ideal ones given in equation (27), 
their effect on increased error may be small. We certainly cannot tell 
from the equations just how sensitive the minimum error will be to changes 
in the fixed network from the optimum form in (27) since, as we have seen, 
we could not even tell whether the effect would be beneficial or detri- 
mental. However, it has been shown in reference 2 that the effect on 
minimum error of the dynamic factors in equation (26) is small as long 
as the natural frequencies are not too low and the damping ratios are 
reasonably small. A good many vehicles fall in this category. 

Solution .- With these forms for \|r s (u)) and Hp(co) , it is possible to 
solve equation (19) exactly. This exact solution is given in appendix A, 
and it can be seen to be quite unwieldy. For practical purposes we would 
like to know if suitable simplifications can be made without sacrificing 
appreciable accuracy. 

Such simplifications can be made as discussed in appendix A. The 
nature of this solution will now be outlined. First of all for the forms 
of i|f s (u) , t n (w) , and Hp(to) given in equations ( 2 l) , ( 22 ), and (27), the 
error equation (l 4 ) , the restricted quantity equation (15), and optimum 
compensating network equation (19) can be made dimensionless by means of 
the following substitutions: 
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CO = Px 

p = %faj^ 

v = s/p 


(28) 


Obviously, the parameters (3 and v are associated entirely with the input 
characteristics (i.e., signal and noise). The dimensionless analog of 
equations ( 1^+) , ( 15) , and (19) become 


W 1 
NP 



^ 4~ Y 2 X ^ 2 ^ + f I Y ( X ) | 2dX 

X 4 (x 2 +V 2 ) i 00 


R 2 k f 2 = A 2 . = f°° lY(x) l 2 
Np 5 WP 5 x 2 +v 2 


Hco(x) 



dx + 


Y(x) | 2 x 4 dx 



H f (z)\|f s (z)e lpzt 
A~( z) 


P dz dt 


(29) 

(30) 

(31) 


A 

2 

7 


where 


A(x) = 


Hf (x) Hf (x) + p 


[ , 'l r s(x) + tn( x ) 3 


Note that the left side of (30) becomes R 2 kf 2 , where k f is the vehicle 
gain. Since the fixed network has only two poles, both at zero, the 
quantity Rkf is merely the acceleration. Thus 


A 2 = 


n 2. 2 
R kf 


(32) 


It should also be noted that the functional form of the quantities in the 
above equations is changed by the transformation (28) . 

The solution of equation (31) for the optimum compensating network 
Hcc/w) is somewhat involved. For this reason the details of the deriva- 
tion are discussed in appendix A. It is shown that if 


!2- « 1 (33) 

27 

which is satisfied for nearly all cases of interest, the approximate opti- 
mum closed- loop transfer function, Y 0 , can be expressed in dimensionless 
form as follows: 


ro > 
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v / v -Px 2 + i \[2. (\[2 +t])x + 1 . . 

Y °( x ) = 7" ; — 777 — 2 77 — 7-7 — 7 = Y 0 (x;vi,ti) 

( ix+l) ( -x 2 +ix+l) ( -T] 2 x 2 +i s[2r\x+l) 


(34) 


or 


Y 0 (s) = 


P P \f2.(s]2+r\) 

— - s + s + 1 

Pi 2 Pi 


+ 1 


+ — + 1 


Pi 2 1 Pi 


Pi 2 


n/2t1 

Pi 


^ S 2 + s + ! 


(35) 


P .2 . ^2(^2+n) 

Pi 2 s + Pi 


s + 1 


. V2-n(l+ nTSti) . (1+n/2ti) 2 (72+n ) 2 . j2(-j2+r\) 

S D H H S° H 


Pi 2 


+ 


S + 1 


(36) 


where 


P = P(vi,T]) (see eq. (A84) ) 
7 = p 1/4 kf -1/2 
n = 701 


(37) 


The corresponding open- loop transfer function p 0 is 
1 -Ex 2 + i 72(-/2+n)x + 1 


Po( x ) = - 


(•/2+r,) 2 - P 


iT ) 2 


3 ^(1+^211) 2 i(l+72tl) 2 

x - — — -7 x + — = — — X + 1 


(v/2+t) ) 2 - P (72+ti ) 2 - P (72+n) 2 - p 


= p 0 ( x J v i^) 


(38) 


Po( s ) = 


(n/2+t)) 2 - P 




3 [ (^/2+tj) 2 - P] 


>i p 72(72+1]) 

"2 S + A 

Pi 

72 ti( 1+ 72 ti) 
Px 2 [ (72+ T,) 2 - P ] 


S + 1 


(1+ 72ri) : 


Pit (72+ T)) 2 - P] 


s + 1 


(39) 

It will be seen that this solution for Y 0 (x) and Uo(x) is dependent on 
only two dimensionless constants Vi and r] . (The subscript 1 has been 
attached to v and (3 in order to associate these quantities with the 
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inputs for which the system was optimized.. Later a subscript 2 will also 
be used to associate these quantities with the actual input which may or 
may not be the same as that for which the system was optimized.) It is 
clear from the definition in (28) that v x is dependent solely on the 
target characteristics. The other parameter, q, is dependent on the prod- 
uct of the input parameter Pi and the vehicle parameter 7. The param- 
eter 7 is clearly associated with the vehicle from the above definitions; 
that is, 7 depends on the vehicle gain kf and on the Lagrangian multi- 
plier p which is used to place the desired restriction on vehicle accel- 
eration. The determination of p is discussed in later paragraphs. 
Obviously p = 7 = q = 0 corresponds to the Wiener case with no restric- 
tions, and in this case the transfer function simplifies considerably. 


Having found the optimum system one can now determine the error and 
acceleration. Because we are interested in the performance of optimum 
systems for a variety of inputs , we must return to the more general 
expressions given in equations (ll) and (12) . 

The error equation (ll) is seen to consist of four terms which are 
defined to be, respectively. 


The first three terms are due to the signal while the last is due to 
noise. It is shown in appendix B that each of these components can be 
expressed in dimensionless form as follows 


Performance Equations of Optimum Systems 


E 2 = E1 2 + E 2 2 + E 3 2 + En 2 


(40) 


( 



00 1 


'| 1-Y 0 (x) + x 2 Q I 2 
x 4 (x 2 +v 2 2 ) 



— - = 4itiQ(a 2 +Qa 4 ) = f2(vi,v 2 ,T)) 


(42) 




2 


= 2itQ 2 = f 3 (vi,v 2 ,q) 


(43) 
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■n 2 f~\ 00 

iitl = / |Y ° (x) |2(Jx °"‘ 1 T t - (44) 

—00 

where Ai, A^, A 3 , and A 4 are determinants given in appendix B and Q 
is a parameter yet to he explained. 

There are several general comments to he made concerning these 
equations. First, it will he seen that all the components of error can 
he readily evaluated hy means of simple algebraic expressions or deter- 
minants. Second, the subscript 1 has been attached to certain quantities 
to identify them with the input for which the system is optimized, while 
the subscript 2 has been attached in order to identify these quantities 
with the actual input. Third, it can be seen that all the dimensionless 
components of error are functions of three parameters, Vi, v 3 , and rj . 
Actually for many cases, these components are functions of only two param- 
eters. For example, for optimization problems (where the input is the 
same as that for which the system is designed) , it is necessary that 
Vi = v 3 * In other problems, where the design input and the actual input 
are not the same, we will see that either v x or v 3 is zero. In these 
cases, the components will be functions of only two parameters. Fourth, 
it should be noted that the total error cannot be obtained in dimensionless 
form by adding equations (4l) through (44) since the nondimensionalizing 
factors are not all the same for all components. And last is the factor 
Q which is somewhat involved and needs some explanation. From the defi- 
nition in equation (B24) , 


Q = f / lim (45) 

S-*0 S ( S+ ^2) 

and (36) for the optimum Y Q , it can he shown that Q will always he zero 
■unless Vi ^ 0 and V2 = 0, and in such a case its value will he dependent 
solely on Vi and tj. That is, 

Q = _a 3 - P = Q( Vi,T]) if V! 0, v 2 = 0 -i 

\ (46) 

= 0 otherwise J 

The physical meaning involved here is simple. For an accelerating target 
at least two integrations are always required in the forward paid: of the 
open-loop system. As can be shown from equation (36) for the optimum Y 0 , 
the least number, namely 2, occurs when ^ 0, or when the system is 
optimized for a signal of the form \|r^(u)) = a/(co 2 +v 2 ) . In such a case Q 
will not be zero only if the signal input to this system is of the form 
cr/u) 2 , that is, = 0. A nonstationary input of this form will be 
discussed in the next section. 

Consider now evaluation of the restricted quantity given in 
equation (l2) . It is indicated in appendix B that when Hf and c r are 
substituted in equation (12) , the vehicle gain kf will combine with R 
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so that the left side of (l2) will become R 2 kf 2 . Since there are no 
poles of the fixed network other than at zero the quantity Rkf is merely 
acceleration. Equation (l2) can then he expressed as 


A 2 = A x 2 + A 2 2 + A 3 2 + A n 2 
= A s 2 + A n 2 


( 47 ) 


where the signal component is 

A s 2 = A x 2 + A 2 2 + A 3 2 


( 48 ) 


Each of these components can he written in dimensionless form as 
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N 2 (3 2 6 


>° °iy(x) - 

x 2 + v 2 2 


, jti ^5 / \ 

ax = - 


( 49 ) 


Ao 2 = 0 


— 


*n 


P s Tx 2 2 % 2 


jti 4 6 


N 2 P: 


x 4 | Y(x) | 2 dx = — 1 r~ = g 2 (v 1 ,ri) 


a o A 4 


( 50 ) 

( 51 ) 

( 52 ) 


where again the V s are determinants and are given hy the complete 
expressions in appendix B. As before, one can show that the components 
are functions of three parameters v x , v 2 , and T| , although for most cases, 
the components will he functions of only two parameters. Also, the total 
acceleration cannot he obtained in dimensionless form hy adding equa- 
tions (49) through ( 52 ) . Later, however, when we have reason to put 
Pi = P2 (or simply P), the components can he added. The parameter M A 
is quite similar to the previous parameter Q. From its definition 


„ . sY(s) n . s 

Ma - lim — = lim 

s-*o ( s+ £s) s-^o s + ^2 


it is clear that 


M A = 0 
= 1 


if v 2 ^ 0 
if v 2 = 0 


( 53 ) 


( 54 ) 


Physically, this condition is obviously dependent on the nature of the 
input signal. 
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Resume 


Because of the length of the previous sections , a resume will now 
he given in order to summarize the problem, assumptions, and solution. 

The class of guidance systems considered above are those described in 
sketch (a) and the accompanying discussion. Our primary objective in this 
report has been discussed in the Introduction. Briefly, it is to relate 
performance, both optimum and off-design performance, explicitly to the 
guidance and control task, that is, to those parameters which are normally 
included in the statement of the problem. 

Since a system cannot be optimized for all inputs and all restrictions, 
it is necessary to specify certain quantities. It is necessary to know 
something of the signal and noise inputs, the number of saturating ele- 
ments, and the fixed network (see eq. ( 19 ) and succeeding discussion) . 
Although a solution can be obtained for any choices of these quantities, 
we would like to choose them so that we can get an explicit solution and 
yet have these choices be as physically meaningful as possible. The 
choices which were made are: 

(1) The signal is of the form given in equation (20) . This form 
fits several inputs, as will be seen in following sections. 

(2) The noise is stationary, white, and uncorrelated with the 
signal (see discussion following eq. (22)). 

( 3 ) Only the most critical saturating quantity (acceleration) 
and one fixed network need be considered (see p. 9 of ref. 1 for jus- 
tification) . Furthermore, even though the fixed network is generally 
very complicated, it can be adequately represented by Hf = kf/s 2 
(see the discussion in connection with eq. ( 27 ), and ref. 2) . 

There are, in addition, several assumptions made in obtaining the explicit 
solution: 

(4) The system response time is less than the interval of flight 
time 0 to T (to be justified in a later section) . 

( 5 ) The weighting factor P is independent of t. 

(6) The following inequality, which is not really essential, 
greatly simplifies the expressions (see eq. ( 33 ) and- discussion) . 
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The solution consists of two parts: the optimum transfer functions, 
and the performance equations for the optimum systems. The dimensionless 
optimum transfer functions are given by equations (34) and ( 38 ) . It is 
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seen that "they have "been expressed, in "berms of "two dimensionless 
parameters, v x and n (defined in eqs. ( 28 ) and (37), and described below 
eg.. ( 39 ) ) . The transfer function (34) is optimum for several inputs and 
forms the basis for optimum time-varying homing systems (as will be shown 
in later sections) . 

The performance equations for the optimum system consist of the 
error equations ( 4 l) through ( 4b) and the acceleration equations (49) 
through (52) . All of them are given in dimensionless form as functions 
of the three parameters Vi, v 2 , and rj. These equations can be used for 
constant-coefficient systems and time-varying homing systems, and for 
certain stationary and nonstationary inputs. They can also be used for 
evaluating off-design performance. Such uses of these equations are 
illustrated in following sections. 


OPTIMUM PERFORMANCE 


A 

2 

7 

4 


It will be the purpose of this section to show how the solution just 
derived applies to several distinct types of optimization problems, and 
to present the corresponding theoretical optimum performance curves. The 
problems considered involve different types of systems and inputs, both 
time-invariant systems and time-varying homing systems, and several types 
of signal characteristics both stationary and nonstationary. 


Time-Invariant Systems 


There are many guidance and control problems that can be described 
by constant— coefficient differential equations with stationary signal . 
inputs which possess frequency characteristics of the form (20) . . It is 
of interest to enumerate some random processess which fall in this cate- 
gory. All have been described previously elsewhere, although in difierent 
terms. When they are expressed in terms of the same definitions, these 
processes and their descriptions can be summarized as in figure 1. It 
will be seen that in signals A and B, both the amplitude and the interval 
length are random variables, while in C only the interval is random. 

The optimum performance curves for systems subjected to signal inputs 
of the form ( 20 ) can be readily obtained in dimensionless form from the 
performance equations presented in the last section. For this case let 
us see what some of the parameters should be. First, since we are con- 
cerned here only with optimization problems, the input to the system is 
to be identical with that for which the system is optimized. In this case 
0! equals p 2 and- v x equals v 2 , and we may therefore use simply P and v 
(also N with no subscript) . Furthermore, since the system is to be 
optimized for inputs of the form ( 20 ) , it is necessary that v not be 
zero. Second, it is clear that the form for the signal input given in 
equation ( 20 ) is applicable to all the processes in figure 1 provided the 
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symbols p and | are properly interpreted in terms of the particular 
process of interest. Third, for stationary processes the input weighting 
function P^ in equation (2) is 2jt; that is, P s = Pjj = 2rt. Fourth, the 
factors Q and % which are involved in the performance equations can he 
seen from equations (46) and (54) to he both zero. To summarize the param- 
eters involved, then, we must have in the performance equations 


Vi = v 2 = v ^ 0 

Pi = P2 = P > 

M A = Q = 0 


( 55 ) 


As a result of these values, several components can he dropped from 
the performance equations since 


E 2 — E3 — A3 — 0 

Also, since Pi = P 2 = P "the components of error and acceleration can he 
combined to give 


E 2 _ 

E1 2 Enf 

( 56 ) 

Np “ 

NP 

Np 

A 2 

A1 2 A n 2 

( 57 ) 

ID 

% 

3 

01 

NP 


Each of these components is still given hy the equations (4l) , (44), (49), 
and (52) . It is only the labeling of the left side of these equations 
which has been changed since P does not occur on the right. Plots of 
each of the components of error and acceleration as well as their total 
are presented as functions of v and p in figures 2(a) through 2(d) . The 
dependent quantities may he taken to represent either ensemble or time 
averages . 

The curves presented in figures 2(a) through 2(d) may he used to 
evaluate optimum performance for any set of conditions for which 
4v s /27 « 1. Given any signal or target parameter and noise as defined 
"by v, P, anfl N, and a vehicle with any rms acceleration capability A, 
the factor t] can he found from figure 2(d) . From this rj, the minimum 
error can he determined from figure 2(h) , and the optimum system transfer 
functions from equations (35) and (39) . If, however, one is not interes- 
ted in knowing the optimum transfer function, the intermediate parameter 
■q can he eliminated as shown in figure 2(e) , and the minimum error can he 
obtained directly as a function of available vehicle acceleration. 
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The impulse responses of the optimum systems are of fundamental 
importance. They are useful in establishing minimum launching ranges; 
they determine the error or miss due to an impulse of signal or noise; 
and they determine the minimum duration of certain nonstationary inputs 
for which the solution is valid, as will be seen. The impulse responses 
are shown in dimensionless form in figure 3* Note that v has a much 
smaller effect on these responses than does q. 


Time-Varying Systems 


In this section is considered a class of important control system 
problems which is characterized by time-varying differential equations . A 

The nature of this class of problems is illustrated in the following sketch 2 
by the example of a homing missile intercepting a bomber target. It is 7 

seen that this example is characterized by the fact that the range between 4 



the missile and target changes continuously with time. This range varia- 
tion is due to physical facts and hence cannot be avoided. Since range 
enters into the coefficients of the differential equation describing this 
problem, the problem is necessarily time varying. The time-varying situa- 
tion illustrated here may be recognized as belonging to a large and impor- 
tant class of problems such as mid-course or terminal guidance in 
interplanetary flight, fire control, and aircraft landing. 

This time-varying problem has the same ingredients as does the 
time-invariant problem just discussed, that is, the target maneuver, the 
noise, and the missile maneuverability. In addition, however, the system 
is constrained to operate with a forced time variation representing the 
varying range. Furthermore, in this type of problem we are concerned with 
minimizing the error only at a particular time T (the time of arrival at 
the destination) , but restricting the vehicle capabilities at all times 
t 2 previous to T. Thus ensemble averages are particularly meaningful. 
From these remarks then one can see that this time-varying problem is 
similar in concept but basically more complicated than the time-invariant 
case discussed in the previous section. 

It is not at all obvious that the solution for the time-invariant 
case given in the previous section has any connection with the present 
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■time-varying problem; however , a relationship can be shown. The 
optimization of the time-varying problem has been studied in references 1; 
a.nrl 6. There an equivalence was shown between the time -invariant and the 
time-varying homing problems, so that the solution given in the present 
report forms the basis for the optimization of the time-varying system. 
More specifically, it is clear from the results of references 4 and 6 
that: 


(1) The optimum performance curves (for the rms ensemble average 
of the error at time T, or miss, and the rms ensemble average of the 
acceleration) given in this report in figures 2(a) through 2(e) are 
valid for the time-varying problem. 

(2) The synthesis of the optimum homing system can be accomplished 
by combining the results of this report and those of reference ^ for 
the homing study. That is, the transfer function Y 0 (s) given in 
equation (35) of this report can be substituted for H(s) in equa- 
tion (22) of reference 4. From this latter equation one may then 
synthesize an optimum control system which, incidentally, is also 

t ime - varying . 

It is not intended to deal further with the details of the optimization 
of the time-varying homing systems since this problem was the subject of 
references k and 6. 


Nonstationary Inputs 


One of the assumptions generally made in the Wiener theory and in 
Newton' s modification of this theory is that the input process must be 
stationary. However, in many physical problems, especially in the inter- 
ception of targets, the input process must be considered as essentially 
nonstationary. The reason is that all real target maneuvers will have a 
finite beginning and a finite end. Consequently, it will be desirable to 
examine the applicability of the solution presented earlier to certain 
nonstationary cases of interest in the interception problem. 

Nonstationary step signal.- One type input of interest in the target 
interception problem is a step of target acceleration as is indicated in 
the following sketch, where the beginning of the step is equally likely 


Target H | h 

acceleration q j 

t — - 


to occur anywhere within the interval of interest. This interval starts 
at some finite time, t = 0, when the vehicle is launched and extends to 
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some time T when the vehicle reaches the target. Such an input is a 
random process which is distinctly nonstationary. It should he noted that 
this type of maneuver has "been used previously in reference 7 • Optimiza- 
tion for this input will he shown to he a special case of the solution 
presented earlier. As will he seen, the solution given here is one which 
enables the results for this and other inputs to he unified. 


Let us consider first the problem of merely evaluating the mean- 
square ensemble average of the error at time T due to a signal only. 
In the real case of interest, of course, a noise signal would he added 
to this input. The real time block diagram with a pure signal input is 
shown in the following sketch where the time ti at which the input 


JL 

U 0 (t - t,) 



impulse occurs is uniformly distributed in 0 to T. This diagram is 
obviously of the same general form that was used previously in sketch (h) . 
We have 


cp (s) = af(s) = 

s 

(58) 

W s (s) = 1 - Y(s) 

(59) 


where Y(s) is the over-all closed- loop transfer function. Now let us 
see what the values of the parameters involved in the solution will he. 
First, using equations (17) , ( 58 ); and the fact that P s = l/T (since 
ti is uniformly distributed) , we have 


*s(“> » t 60 ) 

2rtTu) 

This latter function is certainly not a spectral density and can only he 
called a frequency function which is associated with the signal. However, 
one can see that this factor for the nonstationary process enters the 
equations exactly as would a spectral density for a stationary process. 
Hence one can derive minimum errors and optimum system transfer functions 
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even though the process is distinctly nonstationary. Second, hy comparing 
equations (60) and (2l) we see we need only make the following definitions. 


a 



2*T 


(61) 


i = 0 


Third, since the system is to he optimized for this input, it is necessary 
that £1 = 0, and since the system is to he subjected to this input, 

| 2 = O. In dimensionless terms (see eq. (28)), this amounts to 


Vi = V 2 = V = 0 


(62) 


Fourth, from equations (46) and (54) we see that 


Q = 0 
% “ 1 


(63) 


Fifth, since the actual input is the same as that for which the system 
is optimized, 


Pi = P 2 - P 

The performance curves can now he readily obtained from the solution 
given in equations (4l) through (44) and (49) through (52) . Note that 
again the E 2 and E3 components of error are zero; also A 3 = a 2 = a^ 
and is therefore not a function of v x , v 2 , or T). Thus the performance 
equations simplify to 


E 2 E x 2 En 2 

Np NP + NP 


(64) 


A 2 - a T 2 A x 2 A n 2 

Np 5 Np 5 NP 5 


(65) 


The components on the right are still given hy equations (4l) , (44), (49), 
and (52) , and they are now functions of only one parameter, t), since the 
v' s are zero. 3 At this point it is both illuminating and useful to 
examine the transition that occurs from the stationary case just presented 


3 Note that the Ei and A x components given hy equations (B37) and 
(B72) appear to he indeterminant. However, that this is not the case is 
clear hy expanding Ai, A 2 , and A 5 along the last row. Hence one can: (l) 
simply eliminate the last row and column of Ai, A 2 , and A 5 or (2) merely 
use a very small value of v 2 , say 10" 4 . 


24 


to this nonstationary case. From the details of the equations vhich are 
presented in appendix B one can deduce that Ex, E n , (and therefore E, 
the total error) and An components are part of a continuous transition 
from the stationary case. That is, all of these components will be iden- 
tical with that for the stationary case for Vx = v 2 - 0 (see figs. 2(a) , 
(b) , a nrl (c)). For clarity, plots of these components are repeated in 
figures 4(a) and (b) . The remaining two components, A x and A 3 , comprise 
the total signal component of the acceleration and their values are given 
by the appropriate equations. The transition for these two components is 
not as simple but can be explained as follows. In the stationary case 
when Vi = v 2 ^ 0, the A 3 component was zero. As Vi and v 2 approach 
zero, the other component A x grows without bound because the first 
integral in equation (30) is improper in the limit. Hence it cannot be 
plotted on figure 2(c) . Physically, this is associated with a final 
steady value of the acceleration. In taking this into account as was done 
in the derivation one gets a different equation for A x (in which % = 1 
instead of zero) and in addition another component, A 3 . In other words 
as we progress to the nonstationary case, the A3, for. the stationary case 
becomes, in the limit, the two terms A3, and A 3 . As shown in equa- 
tion (65) the new term A 3 (which is simply a T ) can be moved to the left 
side of the equation since it is not a function of v or q. The A x com- 
ponent is plotted in figure 4(b) . The total error, equation (64) , and the 
total acceleration, equation (65) , a-re plotted in figures 4(a) and (b) . 

If q is eliminated, the data in these two curves can be cross-plotted 
to give figure 4(c), that is, the minimum dimensionless error directly as 
a function of dimensionless acceleration. These curves may be used to 
evaluate optimum performance in precisely the same manner as described on 
page 19 except that they are now much simpler because only the parameter 
q is involved. 


The optimum system transfer functions for this case are still given 
by equations (35) aud (39) if we put Vx — 0. We have 

,2 


Y 0 (s) = 


(%/2+q)^ \/2(\/2+q) 

5 — s 2 + s + 1 

Pi 2 Pi 
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Pi 


S + 1 
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S 2 + 


s/2q 


( 66 ) 


Px(l+^q) 


S + 1 


(6 7 ) 
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.... The impulse responses of this system described by equation (66) are given 
in figuz*6 3 *by "the Vi = — 0 cujrvGSj ■fch.sir significance has "been 
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discussed previously. One can note that the form of the open- loop 
transfer function is now different; that is , there are now three integra- 

• tions required in the forward loop rather than only two. This difference 
is readily explained. The P which was previously very complicated and 
given by equation (A84) approaches a very simple value , P = (*J2+t\) 2 , as 
Vl approaches zero. This has the effect of increasing one of the time 
lags in the forward part of the loop, and in the limit it becomes an 
integration — a first-order lag term with infinite time lag. The physical 
reason for this extra integration is quite simple. It is due to the 
assumption about the input, that is, that there will be one and only one 
switch of the acceleration. Actually, this can never be stated with cer- 
tainty. Saying that the target definitely maneuvers only once during the 
flight interval is certainly not the same as saying it is not likely to 

A maneuver more than once, and this difference leads ultimately to the 

2 difference in the number of integrations required of the optimum guidance 

7 system . 

4 

Nonstationary signals derived from stationary ones .- Other types of 

* nonstationary signals of interest are those derived from stationary sig- 
nals. The situation is illustrated below where the top sketch shows the 
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stationary maneuver which extends in both directions to infinity. The 
peal maneuver, however, will necessarily start at some finite time called 
zero and will end at some finite time T ‘when the vehicle reaches the 
target, as indicated above. This process might be termed stationary in 
the interval 0 to T, since it is part of a stationary process. 

Since the upper of the two inputs does not occur in nature, it is 
often stated that this input is unrealistic and cannot be used. It is 
true that the stationary theory is applicable, strictly speaking, only to 
the upper of the two inputs above. However, the nonstationary character 
of the lower input is due to the mathematical definition of stationarity . 

In the practical case it is clear that it makes little difference to the 
vehicle, so far as error is concerned, whether the process persists over 
an infinite or finite period so long as the process begins at a time before 
the interception point by an amount equal to or greater than the system 
response time. (Of course, the process may terminate any time after time 
T without affecting the results.) In other words, an infinite period is, 
for practical purposes, simply one which is longer than the system response 
time. Thus when response times are short, results presented previously 
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in figure 2 for time-varying and time-invariant systems apply to the 
nonstationary case cited. Fortunately most interception situations fall 
in this category. The impulse responses which were presented in 
generalized form in figure 3 can be readily used to verify this condition. 


OFF-DESIGN PERFORMANCE 


In previous sections we have considered only optimization problems, 
that is, problems in which the system is subjected to precisely the same 
input for which it was specifically designed. Here we will consider off- 
design performance, that is, the deterioration in the error when the system 
is optimized for one input but subjected to a different input. It is 
important to note that the actual input might be different from the design 
input for two reasons. First, the two inputs, the actual and design 
inputs, might be describable by the same type of process, but the numer- 
ical values of the parameters for each input might be different. Such 
situations would occur, for example, if a noise level were different from 
the design value, or if there were a change in signal magnitude. Second, 
the actual and design inputs might be different because the type of 
process is different. 

In the following sections different off -design cases will be 
considered. The manner in which the previous solution can be applied to 
these cases will be immediately apparent. However, to avoid confusion 
between the various parameters and components of error and acceleration, 
each section will be organized thusly: first the parameters involved in 

the solution are given; second, the error and acceleration components are 
arranged in tabular form together with the figure number for components 
which are plotted. It should be noted that the components are given 
individually since they cannot be combined because of different nondimen- 
sionalizing factors. Where possible, the components are combined. It 
should also be noted that certain curves will be identical with previous 
curves. Rather than repeat these curves it will only be necessary to 
alter the ordinate to the dimensionless form indicated. In so doing, the 
ordinates of the curves must be expressed as rms values. 


Actual and Design Inputs of Same Type, Values Different 


Stationary inputs.- Let us assume that a system has been optimized 
for any of the stationary signals (in addition to noise) which have been 
discussed previously. Then, let us assume that any (or all) of the spe- 
cific values of the input process for which the system was optimized are 
changed. The solution already obtained is applicable to these problems. 
Since the actual input and design input are both the same type of station 
ary process, but the numerical values describing the processes are not 
necessarily the same, we must have 
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Vi / V 2 / 0 

Pi 3* p 2 


Also from equations (46) and (54) we see that 


Q = 0 


M A = 0 


The error and acceleration components for this case can now be 
summarized as follows: 



fl( v l^ v 2^ T l) 


not plotted; see equation (B37) 


— 0 
E3 = 0 

= f 4 (vi,T]) rms values plotted in figure 2(a) 

%Pi 



g 1 (vi,v 2 ,n) 


not plotted; see equation (B72) 


A3 — 0 

= g (vi,tj) rms values plotted in figure 2(c) 

N 2 Pi 5 


It can be seen that the E x and A x components are functions of the three 
variables Vi, v 2 , and t \. To display these components with adequate 
accuracy, a good many curves would be required, and they are therefore 
not plotted. (This is the only case in which three variables are required.) 
In a specific application the desired curves could be readily obtained 
from the equations indicated. The other components are functions of only 
the two variables, Vi and rj . 
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Nonstationary inputs.- Let us now see how the solution applies to 
the situation when the system is optimized for specific values of the step 
maneuver and noise , hut the values of the actual signal and noise are dif- 
ferent. In this case we must have 

v-j_ = v 2 = 0 Q = 0 

Pi ^ Pa % = 1 


From these values, the performance curves can be obtained from the 
solution given in equations (4l) through (44) and (49) through (52) . The 
components are 
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e 3 = 0 


E n 2 
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A 3 2 = a T 2 

-^4 = sM 
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rms values plotted in figure 4(a) 


rms values plotted in figure 4(a) 


rms values plotted in figure 4(b) 


rms values plotted in figure 4(b) 


We see that since the A 3 component is a constant, it has been combined 
with the total signal component As as follows 

A s 2 - a T 2 = A3 2 

In this way, total acceleration can be obtained directly from figure 4(b) . 
Also note that all the components (other than A 3 ) are functions of only 
the one variable t] . 
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Actual and Design Inputs of Different Types 


System optimized for stationary signal, actual signal nonstationary . - 
Let us assume that a system has been optimized for any of the stationary 
signals (in addition to noise) which we have discussed previously, hut 
that this system is to be evaluated against a different type of maneuver 
such as the nonstationary single step maneuver also discussed before. The 
previous solution can be shown to be immediately applicable . From the 
discussion of earlier sections it is clear that we must have 


v x ^ 0 Q i 0 

v 2 = 0 Ma = 1 

Pi i P2 


For these values the error and acceleration components can be 
summarized as follows 

2 7TI 2 

E i +Ea x = f 1 (vi,n) + £2(^1, T]) rms values plotted in figure 5(a) 

/ ~ M 2P2 6 X 

Pi 5 


**£) (PxT) 


= f 3 (vi,Tl) 


rms values plotted in figure 5 (t>) 


rr- = f4( v i^n) 


rms values plotted in figure 2(a) 


N 2 P 


A s 2 ~ a T 2 _ Al - g (vi,tj) rms values plotted in figure 5 (c) 

/N 2 p 2 6 ^ j ^ 2 p 2 6 ^ 

A 2 2 

A3 — arp 

A r 2 _ g 2 (vi,T]) rms values plotted in figure 2(c) 

W 2 Pi 5 


It is worth noting that now none of the components are zero and that they 
are all functions of the two parameters v x and t). Also, as in the pre- 
vious section, the A 3 component can be combined with the total signal 
component A s . In a later section a specific example will be considered. 
It is also worth noting that the Ei and A x components appear to be 
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indeterminant . This is not the case, however. Possibly the simplest 
procedure (which is readily shown from the equations) is merely to 
eliminate the last row and column of Ai, A 2 , and A 5 . 

System optimized for nonstationary signal, actual signal stationary . - 
For this situation let us assume that the system has been optimized for 
the nonstationary signal (and noise) discussed previously, that is, a 
single step of acceleration occurring any time during the time of flight. 
The actual signal, however, is to be any of the possible stationary sig- 
nals described before. In this case we see that we must have 


Vi = 0 

Q = 0 

v 2 i- 0 

m a = 0 

Pi t P2 



The error and acceleration components now are 


A 

2 

7 

4 


Ex 2 

/N 2 p 2 6 \ 

VPi 5 J 


fi( v 2 ,t}) 


= 0 

E3 = 0 


En 2 

N 2 pr 




rms values plotted in figure 6(a) 


rms values plotted in figure 4(a) 


Ai 


2 


(Sfi) 


8i(v 2 ,ti) 


rms values plotted in figure 6(b) 


A 3 = 0 
A ^ 

— £ — = g 2 (T)) rms values plotted in figure 4(b) 

MaPi 5 

Note only that the Ei and A x components are functions of v 2 and r\ 
while the noise components are functions of only T). In the next section 
a specific example of this case will be considered. 


31 


EXAMPLE OF EFFECT OF TYPE OF SIGNAL INPUT 


The effect on performance of the type of signal input is of particular 
interest in guidance system design. The reason is that one hardly ever 
knows with certainty what type of process the input will he. Consequently, 
two aspects of the problem are important: 

1. The effect on performance of the type of signal input for which 
the system is optimized, and 

2. The effect on performance of subjecting a system optimized for 
one type of input to an input of different type. 

The two problems posed above really amount to a comparison of inputs. 
The question is: What is a good way to compare two inputs which are dif- 

ferent types of processes? Unfortunately, the answer is not clear cut. 

We can see that since the processes are different and therefore the param- 
eters describing the processes are different, the parameters could be 
chosen arbitrarily and independently. In this case the two effects on 
performance we are examining could be arbitrarily large. The previous 
sections could be used to obtain the exact amount of this difference. 
However, a more reasonable comparison might be based on finding conditions 
under which the inputs are equivalent in terms of performance. It would 
appear that P might be a good parameter on w hich to base the equivalence . 
For one thing, from the definition P = G sja/K it can be seen that p con- 
tains all the information about the input signal and noise. Furthermore 
the dominant modes in the optimum transfer function are determined by p. 
Thus one might expect that different types of inputs which have the same 
P would be approximately equivalent. 

Let us illustrate the latter approach by a specific example. We will 
compare the nonstationary input described on page 21 with the stationary 
random input, Case C. For a realistic comparison, it is reasonable that 
the maximum acceleration should be limited to the same value for both 
inputs. Since the values of p for these two inputs are 


Stationary input C: 



( 68 ) 


Nonstationary input: p 



j 
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we see that these will he the same only if 

l = 4t 

Now let us examine the effect on performance when the above criterion 
is used. It is clear that the results of the previous sections contain 
the desired answers. However , since the nondimensionalizing factors were 
not always the same, it will he necessary to remove these factors hy using 
specific numerical values. Since the equations and curves are dimension- 
less, comparisons could he made for other cases of interest. Let us take 
arbitrarily, the flight time, the maximum acceleration, and the noise mag- 
nitude to he the following values: 


T = 10 sec 


a T = 0-95 g 

N = 15 ft 2 / radian/ sec 


then for the stationary maneuver we would have Z - 40 seconds. For these 
conditions p * 1.0 and Vi — V 2 . - v = 0.0$. 


Now we can answer the first question, the effect of choice of target 
maneuver on which system design is based. This effect is shown by. com- 
paring the performances of two systems each optimized for the two inputs 
just discussed. The performance curves can he readily obtained. from fig- 
ures 2(e) and 4(c) , and the result is shown in figure 7 where minimum 
theoretical error for each input is plotted against the vehicle rms accel* 
eration capability. It can be seen the differences between these two 
curves is quite small over the entire range of vehicle acceleration and 
amounts to only a few feet. Thus the difference in optimizing for these 
two apparently different inputs is small provided the P’s are the same. 
It is of interest to note that this difference in performance is small 
even though the actual change in acceleration for the stationary maneuver 
is twice as severe as for the nonstationafy maneuver. 


The second question, the effect of using an input for which the 
system was not designed, can also be readily answered. By utilizing the 
data in figures 5 and 6, one can readily evaluate the effect of .using the 
nonstationary input with a system optimized for the stationary input and 
vice versa. The result is shown in figure 7 hy the two sets of points, 
rather than curves, in order to avoid confusion. It might be a little 
surprising that the deterioration in error from the optimized curves is 
so slight. 


rO > 


33 


CONCLUDING REMARKS 


It seems desirable to emphasize the viewpoint which has "been taken 
in this report. In general, optimization theory expresses certain types 
of operations to he performed in order to determine an optimum system. 

One would like to he ahle to carry out these operations in order to draw 
general conclusions about the best theoretical performance which can he 
achieved and the design of optimum guidance systems from a knowledge of 
the guidance and control task. It cannot he done without some loss, or 
narrowing down of this task. That is, it is impossible to choose a system 
which is optimum for all inputs and all restrictions. Here the desired 
solution has been obtained: first, by restricting the class of input 

signals but yet restricting it to a very useful and important class which 
includes certain stationary and nonstationary inputs; second, by simplify- 
ing the form of the fixed network or output element as was indicated by 
the results of a previous study; and third, by making approximations in 
the analysis so that a simple explicit solution could be obtained without 
sacrificing significant accuracy. The latter two are actually not very 
restrictive. 

There are many ways in which the solutions presented here might be 
used. First, they might be used to determine the best theoretical per- 
formance which could be achieved for any specific case where the vehicle 
and target characteristics are quantitatively known. The result might 
then be compared to the performance of any other system to indicate possi- 
bilities for improvement. Second, the results might be used in preliminary 
design to evaluate the relative importance of each of the factors which 
affect minimum error. Such evaluations are useful in determining those 
design changes which would be worthwhile in attaining smaller errors. Or 
last, the results might be used to reach conclusions about the effect of 
different input signals, or target motions, as was discussed in the last 
section. 

It is believed that the solution given might also be applied to 
other signal inputs not considered here. An example of such a signal 
which is of practical importance in the interception problem is a target 
maneuver consisting of a single switch of acceleration from a negative 
to a positive value at some random time. 

An important extension is needed to the nonlinear problem. It will 
be recalled that the results presented here are based on rms values of 
the restricted quantity, the acceleration. However, most systems have 
hard limits. To insure that the system remains linear, the rms values 
must be chosen small enough compared to the saturation limits . In most 
cases the answer so obtained is near optimum, that is, the rms values 
required to keep the system essentially linear are still large enough that 
there is not much deterioration in error even from the infinite rms value . 
However, there is an increasing number of guidance systems which must 
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operate under more difficult and adverse circumstances in which the vehicle 
maneuverability becomes small. In these cases the error increases rapidly 
and it is important to utilize the maneuverability of the vehicle in a 
more efficient way, that is, by nonlinear control. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Oct. 24, i960 
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APPENDIX A 

DERIVATION OF OPTIMUM SYSTEM TRANSFER FUNCTION 


It will be the purpose of this appendix to derive an approximate 
solution for the optimum system transfer function Y 0 * As described in 
the text we must first find the optimum compensating network H C o- An 
expression for H co was given in the text in dimensionless form by 
equation (3l) and- it is repeated here. 


H co U) = 


2xA + (x) 


f ” H f (z)fe(z)etf£ 3 dz dt 
^-00 A-(z) 


where 


) (Al) 


A(x) = 


H f (x)H f (x)+p 


^ii(x) 


In these expressions x is a dimensionless angular frequency related to 
(*> by 


w = 0 x- 


(A2) 


P = 



(A3) 


where the parameters a and N are related to the input quantities as 
given in equations (21) and (22). The quantity t i;L (x) in (Al) is. 
defined as the sum of signal and noise frequency functions, respectively. 


^ ii ( x ) = i|r s (x) + V n (x) ( A M 


The quantities A+ and A - are defined as the factors of A with poles 
and zeroes in the upper and lower half planes, respectively. Thus 


A = A+A" 


(A5) 
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The input quantities t s ( x ) 811(1 ^n^ x ^ are > from equations (2l) and (22) 


N 


*s(*) = x 4 (x2 ;t 2 -) 


t n (x) = N 


(A6) 

(AT) 


where 


v = — 

P 


(A8) 


It might appear a little odd that N would appear in equation (A6) for 
the signal frequency function, but this is due to the nondimens ionalizing 
factor in (A3). Note also that some of the poles of ty s ( x ) lie on the 
real axis which the theory in deriving equation (Al) does not permit- 
This problem was discussed in appendix A of reference 1, where the correc 
procedure was indicated. The procedure consisted in modifying equation 
(A6) so as to displace the poles at zero slightly off the real axis, and 
after the final answer was obtained letting the magnitude of this dis- 
placement go to zero. Such a procedure is quite unwieldy. In the 
interests of keeping the expressions as simple as possible, the form in 
(A6) can be used if one is careful to remember at the critical points that 
the poles should actually be slightly off the real axis • The remaining 
function in (Al) is Hf, and as has been discussed in the text, it can 
be well approximated by equation (27)- Thus in terms of dimensionless 
frequency, 


kf 

HfM = - p-T 


(A9) 


Having dispensed with explanations and definitions , we will be concerned 
in the following paragraphs with the solution of (Al). 

Let us find some of the functions needed in equation (Al). Starting 
with %i(x) , we can combine (A4), (a 6), and (AT). 


< CVJ 
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t ±i (x) = ^ s (x) + ^ n (x) 

= w x6+v2x4+1 (A10) 

X 4 (x 2 +V 2 ) 

= N p(x)[-p( -x)] (All) 

x 4 (x+iv)(x-iv) 


If we let x™ represent 






p(x) 


3 

II ( x " x pm ) 


m=i 


= x 3 + h 2 x 2 + + "bo 


(A12) 


where 

hs - "( x pi +x p2 +x p3) 

hi = x pl x p2 + X pi Xp 3 + X p2 X p3 

h 0 “ 

Similarly, 

3 

-p(-x) = n ( x+x pm ) 

m= i 

= x 3 - h 2 x 2 + hxX - h Q 




} (A13) 


J 


(Al4) 


Multiplying, we get 


p(x)[-p(-x)] = x 6 + (2b 1 -t> 2 2 )x 4 + ("bi 2 -2b 0 'b 2 ) x2 " ^o 2 


(A15) 
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Comparing equations (A10) and ( A15 ) , we see that the following relation 
between the b’s must exist 


h - 2b 0 b 2 — ® 


(ai6) 


Also it is clear that 


b Q = i 


(A17) 


These latter relations will he needed later. 
Next we know that because of (A9), 


H f (x)H f (x) 


kff_ 

(3 4 x 4 


(A18) 


Therefore 


Hf(x)H f (x) + p = 


x 4 +(kf 2 / pp 4 ) 
x^ 


(A19) 


= P 


q(x)q( -x) 

X^ 


(A20) 


where > if we let Xq n represent roots in the upper half plane,, the 
polynomial q(x) is 


2 

l( x ) “ II ( x ~ x qn) 
n=i 


= x 2 + C-jX + c Q 


(A21) 
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and 


<l( -x) 


2 

n (x+x qn ) 


n=i 


= X 2 - G -j_X + G o 


(A22) 


with 


Obviously^ 


c i — - ( x qi +x q 2 ) 


c o x qi x q2 


> (A23) 


o 2 = 


P0" 


(A24) 


Now combining equations (A 20) and (All) we have for A(x) in (Al) , 


A(x) 


H f (x)H f (x)+p 


tii(x) 


_ „ l(x)q( -x)p(x) [ -p( -x) ] 

- pN — 

x® (x+iv)(x-iv) 


(A25) 


Splitting this expression according to equation (A5) we have 


A + (x) = pN 


q(x)p(x) 

x 4 (x-iv) 


(A2 6) 


40 


A-(x) 


q.( -x)[— p( -x)] 
x 4 ( x+i v ) 


(A27) 


Here we have to recall the remarks made at the beginning of this appendix 
in order to see how to split the x® in (A25) • 

Now the first integral in equation (Al) may be evaluated by 
utilizing equations (A 6), (A9), and (A27), and the result is 


Ii 


r“ Hf(z)faU)e lgZt 
to A ' (z) 


dz 


(A28) 


k# f 00 e ipZt dz 

(B z 2 ( z-iv)q( -z ) [ -p( -z ) ] 


(A29) 



f (z)dz 


(A30) 


This integral is readily evaluated by considering z to be a complex, 
variable and integrating in the upper half plane . We must recall again 
that the second order pole at the origin would actually have been dis- 
placed slightly above the real axis if we had taken the trouble to . o 
so at the beginning of the analysis. Thus there are two poles within 
the contour , a second order pole at the origin and a simple pole at 
z = iv. It is easy to show after a little algebra that 


3 

Res(0) = lim — [z 2 f(z)] 
z-»o dz 


P iv(b 0 c 1 +b 1 c 0 )+b 0 c 0 

vb 0 c 0 t " V*b 0 2 c 0 2 


vb 0 c 0 


(A31) 
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Res (iv) = lim [(z-iv)f(z)] 
z-*iv 


-i(Bvt 

e 

-v 2 [-p(-iv) ]q(-iv) 


= 7 ae 


-0 vt 


(A 32 ) 


Thus I3. in equation (A30) is 


Ii 


2 iri[Res( 0 )+Res(iv) ] 

H 


= — — v 1 


+ 


p 

vPo c o 


t + 7 2 e 



(A33) 


The second integral in equation (Al) is denoted by I2 an< 3 - is 
merely a Fourier transform of Ii« 


X uu 

I le - ll3xt dt 

2jtirac f -(7i+7a)'boCox 2 + (l7iVb 0 c 0 + x + 1 
"boGo (i(3x) 2 (v+ix) 


(A 3 4) 
(A35) 


2 rtiNkf Iq^x 2 + + 1 

4 >o c o 0 2 x 2 (v+ix) 


(A36) 


Note that in ivaluating equation (A3 4 ) it is necessary to carry along the 
pole displacements discussed earlier in order to have suitable convergence 
factors. In the limit as these displacement terms approach zero, one then 
obtains equation (A3 6) . 

Now from equation (Al) the optimum compensating network can be 
found. 
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H co (x) = 


2«A+(x) 


Substituting equations (A 26 ) and (A 36 ) into (A37) gives 


H co ( x ) = 


kf x 2 ( T al 2 x g +2 Cxi^ai x+ ^ ) 

b 0 c 0 (B 2 p p(x)q(x) 


It can be seen from equation (A12) that 


p(x) - n (x-Xpjj) 


m=i 


3 / x 

bo n V"x — + 1 

m=1 \- x pm 


Similarly, from equation (A2l) 


q(x) = n (x-Xgn) 
n=i 


= c 0 n 


n=i 


x 


■Xqn 


■ + 1 


By utilizing equations (A17) and (A24), (A38) becomes 


Hcq(x) = - 


P 2 x2(T al gx2+2^ a , 1 T al x+l) 
kf 3 


n 


m=1 \-xpm / n=1 \-xqn 


2 / x 


+1 n 


+1 


(A37) - 


(A38) 


(A39) 


(a4o) 


(a4i) 
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When we later determine the values of the roots in the denominator, it 
will he seen that some of these factors combine to give 


, v X 2 ( i 2 x 2 +2 £q ilq iX+1 ) 

C ° U " k f (T 3l x+l)(T 71 2 x 2 +2^ 1 T 71 x + l)(T til 2 x 2 +2^ 1 T Rl x + l) 


(A42) 


where 


•pi 


-x- 


•pi 


m 2 _ 1 

T7I - 


2£ 71 T 71 = 


Xp2Xp3 

_ Apg+ X p 3 
\Xp2Xp3 


Tpi 2 = 


x qi x q 2 




2£m*! 


|ii x lii 


x qi +x qg 
x qi x q 2 


Tea 2 = ~(7i+7g)h 0 c 0 


2^aiTai - ( i7iVb 0 Co + v 


(A43) 


J 


Since , as has been shown in the text, the compensating network Hqq 
and the over-all closed-loop system transfer functions Y 0 are related 
by 


Y 0 (x) = H^x^x) 


(A44) 


we can combine equations (A9), (a 42), and (A44) 


Yq( x ) 


Tqx 2 x 2 + 2£ al T a i x + 1 

( Tp !X+1 )( T 71 2 x 2 +2 ^ 7 iT 71 x+l)(T 1 _ Ll 2 x 2 +2^ | _ ll T|_ Ll x+l ) 


(A45) 
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The expressions in (a 42) for the compensating network and in (A45) 
for the system transfer function are the general forms expressed in terms 
of various roots . We must now relate these roots to the various input 
parameters and the restricted quantity in order to get explicit expressions 
for the transfer functions. Let us look first at Xpm which are the 
roots of 

p(x)[-p(-x)] = x 6 + v 2 x 4 +1=0 (A46) 


The substitution 


x 


2 = 


1 

y 


(a47) 


reduces equation (a 46) to 


y3 + V 2y +i=o 


(a48) 


The standard form for the transformed cubic is 


y 3 + py + q = o 


(A49) 


where we see that 


p = v 2 
q = 1 


(A50) 


The general cubic in equation (a 49) can be solved by Cardan’s method- 
However, instead of assuming a solution of the form Y = u + v, it is 
neater to use 


y = v 


u 


ro > 


■F'-nI IU p> 
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Then following the usual development we can show that 

u 3 = A 
V 3 = B 


where 



(A51) 


B 



1 

— + 

2 


2 s/ 27 


(A52) 


and we agree to take the positive square root. Then the first root y 1 
is 


yi = 3 Jb - 3 Ja 


(A53) 


where we "bake "the real roots . Since A is larger “than B, define C as 


c 2 s -y x = 3 Ja - 3 Jb 


(A 54 ) 


Obviously 


-1 < y-L <0 

0 < C < 1 


> (A55) 
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Then we see that 



1 


(A56) 


From (A4-3) 



(A57) 


It is well known that the other two roots, y 2 and. y 3 , of the cubic 
can be expressed in terms of the first root y x . However, it is clear 
from equation (A4-3) that we do not need these roots. It is only neces- 
sary to know Xp 2 Xp 3 and Xp 2 + Xp3 • We may find these quantities from 
(A13) and (Al6); that is, from the four equations 


the four -unknowns, b 2 , b 1 , Xp 2 Xp 3 , and Xp 2 + Xp 3 , can be found. The 
quantities b Q and x pi are known and are given by (A17) and (A56), 
respectively. First, from the last equation in (A58) we get 


b 1 2 -2b 0 b 2 = 0 


- ( X pl +Xp 2 +X p3 ) 


> (A 58 ) 


b x Xp^Xp 2 + XpjXp 3 + Xp 2 Xp 3 


^O _x pl x p2 x p3 



(A59) 


Hence 



1 

C 


(A60) 


ro > 


Now we are left with the first three equations with the three unknowns 
b 2 , "bij and x p 2 + x p3* Solving, we get 


h 2 - 


1 _ 

iC 




(A6D 

(A 62 ) 

(A 63 ) 


where we must take the positive square root in order to get a positive 
damping ratio in the optimum transfer function. From equations (A59) 
and (A 63 ) we can now obtain 


2 ^7i T 7i 


X P2 +X P3 = ± /2 _ — 

X P2 X P3 ^ C 


(A64) 


Next let us find the roots x^ n . From (A19) and (A20), 


q.( x )q( -x) 



(A65) 


Let us define 



ri = 7(3 


(A66) 


Then ( A 65 ) is 


q(x)q(-x) = x 4 + — 

T] 4 


(A 67 ) 
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The roots are 


x qi " k eXp ^?) = J2r) + hr\ 


x q 2 = - « ex P 




n/~2tj + S 2r) 


(a68) 


he needed later. 


we can obtain and c G 

which will 

-i & 

(a 69) 



_i 

(A70) 

Tl 2 


: dynamic factors 


•= • T) 2 

(A71) 

$ 

•H 

II 

(A72) 


The dynamic terms in the numerator of (A42) or (A45) are harder to 
find. Their values are given in equation (A43 ) , and they depend on 7i 
and 7 2 . From equations (A31) and (A32) these parameters are 


7 1 = " 


ivChoCj^+hjCQ) + h Q c 0 

v 2 b 0 c o 2 


(A73) 


(A74) 

- 2 [-p(-iv)][q(-iv)] 


72 = “ 


Let us find 


first . 


Utilizing (A73) we get 
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2£aiT ( 


ai 


2^aiT a i 



(A75) 


Putting in the values of the b's and c's from equations (A 69 ) , (A70 ) , 
(A62) , and (A17) we arrive at 


A 

2 

7 

4 



(A76) 


Note that C is a function of v, so that 2^ iT a i is a function of 
only the dimensionless parameters v and tj. 

The other dynamic term in the numerator of (A42) or (A45) is T ai 2 
given by 


Tq,! 2 _ “(7 1+72 )^o c o 


(A77) 


Putting (A73) and (A74) in (A77), and utilizing equation (A76) just 
developed, we have 


h 0 c o 

T(X1 v 2 [-p(-iv)][q(-iv)] 


iv(b 0 ei+ c obi) 

v2 ^>o c o 


(A 78 ) 


1 _ 

V 2 


"-p(-iv)" 

'q(-iv)" 

L h 0 J 

_ Co _ 


+ i v ( 2 ^q^Tqu ) + 


(A79) 
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Note that all the factors in this expression can be shown to be functions 
of only v and t). This form could be used for computational purposes. 
However, there are two reasons for modifying (A78) further. First, the 
explicit function of v and r| is not displayed. Second, as v becomes 
small, T al 2 apparently becomes large; actually, it does not, since by 
performing some of the operations in the brackets we can show that the 
v 2 cancels out. From (A67), 

-J— + —5^— q(x) 

q(-x) ri 4 x 4 + 1 


Then by use of (A2l) 


1 

q(-iv) 


T] 4 X 4 + 1 


+ 1C ■> V 


+ Co) 


7 

(A80) 4 


Similarly, from (A46) it can be shown that 

p( iv) [-p(-iv) ] = 1 


so that 


— - = p( iv) = -iv 2 — bgv 2 + ibj_v + bQ (A8l) 

-p(-iv) 


Then equation (A78) becomes 




2 


T] 4 CQ 2 ’b 0 2 

v 2 (t| 4 v 4 +l) 




iv 3 


— 

bo 



v + 1 



(A82) 


As shown by equations (A69) and (A70) the c's are functions of T| 
while from (A6l) and (A62), the b's are functions of v. Putting in 
these relations and multiplying, we can eventually show that 


ro > 
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= 


-p 


-1- { L 2 + 

n 4 V 4+l \ 1_ ' 


(c+ -ya 1 ) y § - c 2 + ^T)C + ^ - n 2 c + (n 2 + -/^ic) J § - c 2 + ^f 1 + i 


-n 4 + ^ ^ + n 2 cj| -c 2 v 2 - n 2 - ^n 5 - n 4 (c + /§ - c 2 )jv 3 j (A83) 


Note that only v and ri are involved; C is a function of v as given 
previously "by equation ( A^h - ) . 

For convenience let us summarize all the exact values for the 
dynamic terms of the optimum transfer function which we have just found. 


S*i 2 = 


T (B: 


= iC 


T/i 2 = - 


2£ 7 iT 7 i = i/r- C 2 
V 2 = -T)2 

2^iV = ±^2r 


> 


t] 4 v 4 +1 


n 2 + 


(c+ sf2n) - c 2 + c + ^ 


+ (n 2 + '/2 nc) Jf - c 2 + ^ + 1 


v +n/2ti + £ + n 2 cj§ - c 2 ]v 2 - [n 2 - ^n 5 - n 4 (c +/§ - c 2 ^ 


2 W T ai - 1 


*/2ri + C + ./p - c 2 


(A8U) 
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where 

C 2 = 3 >/a - 3 4b 0 < C < 1 



It will he seen that even though all of the dynamic terms in the 
optimum transfer function are dependent on only two parameters, t] and v, 
the transfer function is quite unwieldy, and should he simplified if 
possible. For this purpose let us look at equations (A51) and (A52) . 

For all cases of practical interest it can he shown that 


so that 


bv 6 

27 


« 1 


A ~ 1 
B * 0 
C * 1 


A 




(A85) 


(A86) 


The dynamic terms can then he shown, to reduce to 


^i 2 = -P 


T pi'= 1 

t 71 2 = -1 

2 ^7i T 7i = 1 

V 2 - -n 2 

2^i_i. = i \[2-T\ 


"N 




— {(n/2+t)) 2 - (l + -J2r\) z v + [v/2t)(1+ J2i\) - r) 4 ]v 2 + [s/2r\ 4 (-J2 +t\) - rj 2 ]v 3 } 

t} 4 v 4 +1 1 


(A87) 


= i V2 (n/2+ti) 


ro > 
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We may now write the approximate closed- loop dimensionless transfer 
function, Y 0 (x ) , as follows. 


-Px 2 +i n /2( n /2+ti)x+1 

Y ° ^ ( ix+l) ( -x 2 +ix+l) ( -r] 2 x 2 +i n/2t[x+1) 


(A88) 


7 

4 


where the P in the numerator is the only complicated term. From (A88) , 
using s = iw = ipx, we have 

x S 2 + {Mhzls + 1 


(A8 9 ) 


n/2('/2+ti) 

t>~ T 

P 2 P 


_P_ g 2 + — - L L s + 1 


+ 


n/2ti(i Win) (lWiriV 


S3 + 


(J2+T]) 2 o2 n/2^1±Ii1 

■ — - — S > ~ " 


s + 1 


p 


(A 90 ) 


N(s) 

D(s) 


One may readily obtain the corresponding open- loop transfer function as 
-Che tLt by equations (38) and ( 3 9 ) • H°te that the dimension- 
less form of the transfer function is a function of only two variables, 
v which is dependent on the input, and n(- P >7 ~ Pp ' which 1 
dependent on the input and the restricted quantity. 
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APPENDIX B 

PERFORMANCE EQUATIONS OF OPTIMUM SYSTEMS 


The optimum transfer function has been derived in appendix A and 
we would like now to derive the performance equations for such an optimum 
system. As indicated in the text we must return to equations (ll) and 
(12). Let us look at the error equation (ll) first. The error consists 
of four components 

E 2 = Ej 2 + E 2 2 + E s 2 + E n 2 (Bl) 


E 


2 

1 



cp s (co)[l-Y 0 (co)] 



(B2) 



2P s c e lim -jq> s (s)[l-Y 0 (s)] 

S-*D 


P s c 6 2 T 


% r 

2 *./ 


c Pn( w )^o( w ) 





(B3) 

(B4) 

(B5) 


Before proceeding to expand the above equations, we will indicate 
certain things done in succeeding paragraphs that are common to all of 
these equations. The optimum transfer function Y 0 (s) in the above equa- 
tions has been shown to be a function of only three parameters: the mis- 
sile parameter q, and the two input parameters and v x . The subscript 
1 is used to associate them with the input for which the system is opti- 
mized. In contrast, all of the parameters other than those in Y 0 in the 
above equations are clearly associated with the input to which the system 
is subjected, and this may or may not be the same as that for which the 
system was optimized. In order to distinguish between these two situa- 
tions, a subscript 2 will be attached to certain parameters to associate 
them with the actual input. Thus equations (20) through (25) would have 
the subscript 2 attached to cr, £, N, and a. Also from (28) 



When the equations are made dimensionless by setting co = p 2 x, as in 
appendix A, it will be seen (e.g., see eq. (B23)) that in equation ( 28 ) 
we would have 



(B7) 
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In the case of 3 X and P 2 the distinction is not very important because 
it will he seen in the text that in many cases of interest we will need 
to put 3i = anyway. The distinction between v x and v 2 is more 
important . 


Let us look first at the noise component E-q. From the discussion 
following equation (2) we know that P n = 2it. Next, we may substitute 
equation (24) into (B5) (using the subscript 2). And finally we can make 
the angular frequency w dimensionless by using the same factor as was 
used in appendix A; that is, by letting w = P^x. Then equation (B5) 
becomes 


Jnf 

N 2 3i 



I Yq(x) | 2 dx 


(b8) 


r in(x)i 2 

x D(x)D(-x) 


dx 


(B9) 


where 

Yo(x) = (BIO) 

Such an expression occurs very frequently and can be evaluated by the 
method of reference 8; that is, an integral of the form 


In 


r g (x) 
SrtiJ^ h(x)h(-x) 


dx 


(Bll) 


where 

g(x) = b 0 x 2n-2 + biX 211-4 + . . . bn-! 
h(x) = a 0 x n + a x x n- 1 + . . . a n 
can be evaluated by the following very simple but elegant expression 


(B12) 


= (-D n+1 % 

2a-o I’n 


(B13) 


where Dn is the determinant given by 


Dn - I ^mrl 

d-mr = a P m-r 0 > 2m-r < n 


= 0 


otherwise 
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and N n is the same as Dn except that the first column of N n is 
replaced hy b 0 , b^ . . .jb n _;L. The N(x) and D(x) in equation (B9) 
can he obtained from equation (3*0 and they are 


N(x) = -Px 2 + i n/ 2( */2+rj)x + 1 
D(x) = ir| 2 x 5 + n/2t)(1+ \[2t])x 4 - i(l+ «/2 t)) 2 x 3 
-[ ( *72+ti) 2 ]x 2 + i-/2( «/2+r|)x + 1 
In evaluating (B9) we have, on comparing (B9) and (Bll) , 

i N(x) I 2 = g(x) 

D(x) = h(x) 

and we see that n = 5. Thus h(x) and g(x) are 

g(x) = b 0 x 8 + bjX 6 + b 2 x 4 + b^ 2 + b 4 
h(x) = a Q x 5 + a x x 4 + a 2 xS + a^ + a 4 x + a 5 
Then it is clear that 

a Q = if] 2 

a. x = \T2t)(1+ •J2x\) 
a 2 = -i(l+ n/2t}) 2 
a 3 = -[( \I~2+t \ ) 2 ] 

8,4 = i\I~2( -f 2+r\) 
as = 1 


(BlU) 


(B15) 


(Bl6) 




(B17) 


J 


b 0 = 0 

bi = 0 

b 2 = P 2 

b 3 = 2( sf2+ r\) 2 - 2P = 2(-a 3 -P) 




(B18) 


J 


A 

2 

7 

k 


b 4 = 1 


I 


The solution of (B9) according to equation (B13), then, is 


E n 
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= = *4^1, T]) 

A 4 


(B19) 


where A 3 and A4 are the following determinants: 


A 

2 

7 

4 


A3 = 


0 

a l 

a Q 

0 


a i 

a o 

0 

0 

b 2 

a 3 

a 2 

a i 

II 

nH 

r< 

a 3 

a 2 

a i 

a o 

b 3 

1 

a4 

a 3 


1 

a4 

a 3 

a 2 

1 

0 

0 

1 


0 

0 

• 

1 

a4 


(B20) 


This is the result given in equation (44) . Since A3 and A4 are simpler 
in determinant form, they will not be expanded. 

Next we will consider the first component of error Ex. From equa- 
tions (8), (10), and (23), 


c e = lim scp s (s)W s (s) 
s->o 


lim 

s-*o 


sa 2 [l-Y(s) ] 

s 2 (s+| 2 ) 


, . l-Y(s) 

= a 2 lim — 15 — ( = a 2 Mg 

s^o s(s+£ 2 ) 


(B21) 


Putting (23) and (B21) in (B2) gives 


E, 2 


P s q 2 2 J' °° I I-Yq(co) + u^Mg - i coggMg 


2it 


w 4 (w 2 +| 2 2 ) 




(B22) 


This latter equation can be made dimensionless, as before, if we let 
oj = PxX. Then 


Ei 2 = P s°2 2 1 


2 jt pn 5 J 

^ -L — O 


I 1-Y 0 (x) + Pi - iP 1 2 M € v 2 x 
x 4 (x 2 +v 2 2 ) 


dx (B23) 


] 
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Let us define Q: 


Q = = Pi 2 lim 

S-K) 


l-Y(s) 

s(s+| 2 ) 


(B24) 


Now if, in the optimum transfer function given in (A45), we put in the 
exact values of the parameters given in (A84) , it can be shown that if 
1^0 (or Vl ^ 0), 1 - Y(s) will have two factorable s's in the 
numerator, while if £ x = 0, it will have three. Thus it is clear that 

Q / 0 if ?i^ 0,| 2 =0 

= 0 otherwise (B25) 

The value of Q if not zero can be deduced from (B24) and ( 36 ) for the 
approximate optimum transfer function to be 

Q = ( -/2+ti ) 2 - P 

= -a 3 - P (B26) 

The last term in the numerator of (B23) is obviously zero. Now from the 

definition of p in equations ( 28 ) and ( 25 ) an alternate expression for 

P 2 is 

P 2 = 6 v/ P s a 2 2 /2itN 2 (B27) 


Then (B23) becomes 

El 2 = f°° ll-Yp(x) + Qx 2 !f ^ (b28) 

^N 2 p 2 6 ^ J -00 X 4 (x 2 +V 2 2 ) 

This latter expression is given in the text as equation (4l) . Putting 
(BIO) in (B28), 

E x 2 _ P I P(x) - N(x) + x 2 QD(x) I ^ (B29) 

/NgPgfN '-[ O0 D(x)D(-x)x 4 (x2+v 2 2 ) 

V Pi 5 / 


By obtaining the N(x) and D(x) from equations (Bl4) and (B15) we can show 
D(x) - N(x) + x 2 QD(x) = x 2 -jiQTi 2 x 5 + Q \/"2t]( 1+ */2q)x 4 + [iri 2 -iQ(l+ n/2ti) 2 ]x 3 

+ [ nT2t]( 1+ - Q( sf2+Tl) 2 ]x 2 + [-i(l+ nT2t)) 2 

+ iQ «J2( ^2+n) ]x - [ ( «/2+t ]) 2 - p - Q]} (B30) 


<$ OJ 
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A 

2 

7 

k 


Thus in comparing (B29) to (Bll) we must have 

h(x) = D(x)(ix+v 2 ) 


g(x) 


D(x) - N(x) + x 2 QD(x) 


x c 


(B31) 

(B32) 


and we see that since D(x) is fifth order , n = 6 . The h(x) and g(x) 
in equation (B12) will now he of the form 

h(x) = eox 6 + e^x 5 + e 2 x 4 + e^x 3 + e 4 x 2 + esx + e 6 (B33) 

g(x) = r 0 x 10 + r^x 8 + r 2 x 6 + r 3 X 4 + r 4 x 2 + rs (B3^) 

After considerable algebra we find that the e's and r's can be 
expressed in terms of the a's just defined in equation (B17) as 
follows : 

e Q = iQ-o 

e i = ia i + v 2 a O 
02 = 10*2 ^2 a l 

e 3 = ia 3 + v 2 a 2 \ (B35) 

e 4 = ia 4 + v 2 a 3 
e 5 = ias + v 2 a 4 
e 6 = v 2 

J 

r 0 = -Q 2 a 0 2 
r x = -2Qa 0 2 

r 2 = -a Q 2 + Q 2 (-2a G a 4 + 2a x a 3 + a 3 2 ) 

r 3 = 2Q[ -a 0 a 4 - a 2 2 + a x a 3 ] > (B 36 ) 

+ Q 2 (— 2a 2 a 4 + a 3 2 ) 

r 4 = 2(a x + Qa 3 ) (a 3 + P + Q) - (&2 + Qa 4 ) 2 
r 5 = (a 3 + P + Q ) 2 

J 


Now we have 
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N 2 P 2 6 ' 
Pi 5 


= 2«il 6 


- - is S - 


where and A 2 are the following determinants 


(B37) 


A-, = 


r o 

e o 

0 

0 

0 

0 


e l 

e o 

0 

0 

0 

0 


e 2 

e l 

e o 

0 

0 


e 3 


e l 

e o 

0 

0 

r 2 

e 4 

e 3 

e 2 

e i 

e o 

a 2 - 

es 

e 4 

e 3 

e 2 

e l 

e o 

1-3 

e 6 


e 4 

e 3 

e 2 


0 

e 6 

es 

e 4 

e 3 

e 2 

^4 

0 

0 

e 6 

es 

e 4 


0 

0 

0 

e 6 

es 

e 4 

1*5 

0 

0 

0 

0 

e 6 


0 

0 

0 

0 

0 

e 6 


(B38) 


7 

4 


Now for the E 2 component of error given by equation (B3) • Using 
equations (23) and (B2l) , and making them dimensionless as before, we 
can show that 

- *}-- # s lvg< * (B39) 

Since the product v 2 Q will always be zero, the last term may be dropped 
and equation (B3) can be written 


E 2 2 = lim X - Y (^ + - ^ - 


-P: 5 


(B40) 


i x-*o x 2 (ix+v 2 ) 

Utilizing equations (BIO) , (Bl4) , and (B15), we get the expansion 

2P s a 2 2 Q ( )x 5 + . . . + [ - i ( 1+ n/~2tj ) 2 + i n/~2 ( n/~2+t|)Q]x - [(n/~2+ri) g - P - Q] 

(B4l) 


e 2 2 = 


-Pi 5 


lim 

x-kd 


D(x) (ix+v 2 ) 


Nov ve can see that if Q ^ 0 (£i / 0 , £2 - 0 from eq. (B25)), the 
value of (B4l) would be 


2 _ 2P s q 2 2 Q 


E 2 2 = 


Pi 5 


(l+ n/2t]) 2 - Q n/2( *7~2+ti ) 


(B42) 
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On the other hand if Q = 0 we can show from equation (B4l) that E 2 = 0. 
For if Q = 0 and v 2 ^ 0, the value of lim | J- will be either 0 or a 
constant , depending on whether V]_ = 0 or not. If v 2 = 0, the denomi- 
nator approaches zero with x; however, according to equation (B25) it 
must be that v x = 0 so that the numerator also approaches zero with x; 
the ratio is a constant. Thus we see that (b 42) is valid whether Q is 
zero or not. Rewriting (b 42) in terms of the a’s used previously we have 

E 2 2 = gPsa a 2Q ( ia 2 + iQa 4 ) (b43) 

3i 


or using (B27) 


E 

2 = 4jtQ( ia 2 + iQa 4 ) = f 2 (v 1 ,v 2 ,ri) 

W 2 (B 2 e 


3 1 5 


which is the same as given by equation (42) in the text. 


(B44) 


The E 3 component is quite simple and can be written in terms of 
the same quantities as the other components. By use of (B2l) and (B27), 
(b 4) becomes 

e = 2 = 


G N 2 T2tcQ 2 

rt A. 


Therefore^ 


e 3 2 


(3iT) 


2jtQ 2 = f 3 (vi,v 2 ,-n) 


(B45) 


All of the error components are now in a dimensionless form and 
are functions of three parameters, v^, v 2 , and r|. As shown in the text, 
for most cases these components reduce to functions of only one or two 
parameters . 

Now we will consider the equation for the restricted quantity given 
earlier by equation (12). The restricted quantity consists of four parts 


R 2 = R x 2 + R 2 2 + R 3 2 + R n 2 


(B k6) 
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where 


—00 

cp es (to)Y(u) c r 2 

sv ' r dto 

Hf(w) iw 

(b47) 

P _ , . rcp s (B)Y(s) C r " 

r 2 2 _ 2P s c r lun H _/ s \ T 

S-*0 L J 

(b48) 

Rs 2 = Ps c r 2T 


(b49) 

R 2 _ ?n r 

Rn _ 2* J 

cp n (w)Y(to) 2 ^ 

■w 

(B50) 


*00 


It is worthwhile to note at this point a common feature of all these 
equations . When Hp and. c-^. are substituted, in the ahove equations } all 
the left sides will "become of the form R 2 kf 2 , where kf is the vehicle 
gain. Since there are no poles of the fixed network other than at zero , 
the quantity Rkf is merely acceleration. 

A 2 = R 2 k f 2 (B51) 


Thus rather than (b 46) we will he interested in the equation 


A 2 = A x 2 + A 2 2 + A 3 2 + A n 2 


(B52) 


Starting first with the noise component , R n , we see that in a 
manner similar to the noise component of error we can express equation 
(B50) as 


Rn 2 kf 2 

N 2 Pi 5 


A 2 

N 2 Pi 5 



x 4 |Y 0 (x) | 2 dx 


(B53) 


Here we have put P^ = 2 jt , used equations (24) and (27 ) } and nondimen 
sionalized hy setting to = IBjX. Then using (BIO) , 


A n 2 r x 4 |H(x)| 2 

(B54) 

N 2 Bi 5 ' io D ( X ) D (’ X ) 

By comparing (B54) with (Bll) we see we must let 


g(x) = x 4 | N(x) | 2 

(B55) 

h(x) = D(x) 

(B56) 


L 
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We have from (Bl4 ) , 

x 4 | N(x) I 2 = -Px 8 + [2( nT2+t ]) 2 - 2P]x® + x 4 


(B57) 


and D(x) is given by (B15). Thus we see n = 5. The h(x) and g(x) 
will be of the form 


h(x) = a 0 x 5 + aj_x 4 + a 2 x 3 + a 0 X 2 + a 4 x + as 
g(x) = iIqX 8 + n x x 6 + n 2 x 4 + n^x 2 + n 4 
Thus the a's are given by (B17) and 

no = P 2 = t > 2 

n ± = 2( */ 2 + t ]) 2 - 2P = b 3 

na = 1 

n 3 = 0 

U4 = 0 


Then we have 


A 


n 


N 2 Pi 5 


= 2trtil5 


- s; r t ■ g2<Vl ’ Tl) 


where 



b 2 

a 0 

0 

0 


a l 

a© 

0 

0 


b 3 

a 2 

a l 

a o 


a 3 

a 2 

a l 

a o 

^6 - 

1 

a 4 

a 3 

a 2 

>• 

II 

1 

9-4 

a 3 

a 2 


0 

0 

1 

9-4 


0 

0 

1 

9^4 


(B58) 


(B59) 


(b6o) 


(b6i) 


Consider next the first component A x given by equation (b 47). 
In a manner similar to deriving equation (B2l) we find that 


c r = lim scp s (s)W s (s) 

S-K3 

* 2 s* (b 62 ) 

kf 
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where 


M, = Urn ^ = Um -V 

s-o s+ 5s s-o s+ ^ 

Obviously there are only two possible values for M^: 

M A = 0 for | 2 ^ 0 

Ma = 1 for ^2=0 

Putting equations (23), (27), and (b 62) in (b 47), and then setting 
cj = 3 x x, we have 


(B63) 


(b64) 


R x 2 k f 2 = A x 2 = 


p s a 2 2 r 

Y(x) 

m a 

2 * lo 

P x (ix+v 2 ) 

i(3 x x 


3 x dx 


Pso-g 2 1_ 
2it (B x , 


ixY(x) - ixMA - M A V 2 I 
x 2 (x 2 +v 2 2 ) 


dx 


P s q 2 2 _1 f 00 lY(x) - m a 1 2 

2jc 3 x J x2+v 2 2 
—00 


(b65) 


Note that Mav 2 is zero in the equation above (B65) • Using (B27) now, 
(b 65 ) becomes 


Ri 2 kf 2 = A x 2 



I Y(x) - MaI 2 
x 2 +v 2 2 


dx 


(b66) 


which is the same as given in the text as equation (49). Now using 
equation (BIO) we have 


A l 2 _ r°° I N(x) - MaD(x) I 2 ^ (b67) 

^N 2 P 2 6 ^ -oo [D(x) (ix+v 2 ) ][D(-x) (-ix+v 2 ) ] 

Using (Bl4) and (B15) it can be shown after a little algebra that 
|N(x) - I^D(x) | 2 = M a 2 ti 4 x 10 + |m a 2 (1+ «/2t)) 4 - M A 2 -/2ri(l+ nTSti) [%( '/2+ 1 l) 2 “ j* x 6 

+ [ M a ( •J~2+r \ ) 2 - P] 2 x 4 + J^2P(M A -l) + 2 (M a -1) 2 (nT2 +1 i) 2 ] x 2 + (M A -l) 2 

(b68) 
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The D(x)(ix+v 2 ) is the same as was used in equation (B3l) , and n = 6. 
Thus in (B12 ) , 


h(x) = e Q x 6 + e^pc 5 + e 2 x 4 + egX 3 + e 4 x 2 + e 5 x + e 6 
s( x ) = ( i 0 x ' 1D + 9i xS + + 93x 4 + d4 x2 + q.5 

Since 


(B 69 ) 


h(x) = D(x)(ix+v 2 ) 
g(x) = | N(x) - M a D(x) I 2 


(B70) 


we see that the e's are the same as given in equation (B35) • The q's 
which are obtained from (b68) can be expressed in terms of the previously- 
determined a's and the result is 


do = -M A 2 a 0 2 
9i. = 0 

q 2 = -M A 2 a 2 2 + 2M A a 1 (M A a 3 +P) 

9.3 = (M A a 3 +P) 2 

q 4 = 2P(M A -1) - 2a 3 (M A -l) 2 

95 = (M A -1) 2 

Finally we have, utilizing the method of equation (B13), 


■\ 


(B71) 


J 



2jtil s 


" §3 T 2 = 


(B72) 
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where 



e o 

0 

0 

0 

0 

0 

e 2 

ei 

e o 

0 

0 


e 4 

e 3 

e 2 

e l 

e o 

0.3 

e 6 

es 

e 4 

e 3 

e 2 

q.4 

0 

0 

e 6 

05 

e 4 

q.5 

0 

0 

0 

0 

e 6 


e l 

e o 

0 

0 

0 

0 

e 3 

e 2 

e i 

e o 

0 

0 

05 

e 4 

e 3 

e 2 

e x 

e o 

0 

e 6 

es 

e 4 

e 3 

e 2 

0 

0 

0 

e 6 

es 

e 4 

0 

0 

0 

0 

0 

e 6 


(B73) 


The second component, A 2 , can he shown to he identically zero. 
Putting (23), (27), and (b 62) in (b 48) we can show 


R2 2 lCf a = A 2 2 = 2Pga 2 2 Mft Hal 

s-*o 


sY(s) - Mfts - %j 2 
s(s+| 2 ) 


Pi 


lim 

x^o 


' Y(x) - % 
iX+V2 


-2P s a a 2 Mft li m |(MA-l) + [iNT2( n/~2+t]) (M A -l) ]x 
Pi lx+v 2 


[Ma(-/2+t)) 2 - P]x 2 + . . .( )xs 




(B74) 


From the relation in (b 6U) we can see that in equation (B7^) 

lim | 1= - i for M A = 0 

x-*o ' ' 

= 0 for M A ^ 0 

Combining (B75) and (B74) it is clear that A 2 2 =0. 


(B75) 


It hardly matters what is done with the last term, A3, since it is 
not a function of the parameters v x , v 2 , or t). Consequently it can he 
moved to the left side of equation (B52) . For uniformity let us sub- 
- stitute (B62) in (b 49) to give 


(B76) 


A 

2 

7 

k 




R 3 2 k f 2 = A 3 2 = P s a 2 2 M A 2 T 


6 7 


REFERENCES 


1. Stewart, El wood C.: Application of Statistical Theory to Beam-Rider 

Guidance in the Presence of Noise. II - Modified Wiener Filter 
Theory. NACA TN 4278, 1958 . 

2. Stewart, Elwood C., Druding, Frank, and Nishiura, Togo: The Effects 

of Target and Missile Characteristics on Theoretical Minimum Miss 
Distance for a Beam-Rider Guidance System in the Presence of Noise. 

NASA MEMO 2-12-59A, 1959- 

3- Laning, J. Halcombe, Jr., and Battin, Richard H. : Random Processes 

in Automatic Control. McGraw-Hill Book Co., Inc., 1956. 

4. Stewart, Elwood C., and Smith, Gerald L. : The Synthesis of Optimum 

Homing Missile Guidance Systems With Statistical Inputs. NASA 
MEMO 2-13-59A, 1959. 

5 . Newton, George C., Jr.: Compensation of Feedback Control Systems 

Subject to Saturation. Parts I and II. Jour. Franklin Inst., 
vol. 254, no. 4, Oct. 1952, pp. 281 - 296 ; and no. 5* Nov. 1952, 

p P . 391 - 413 . 

6 . Stewart, Elwood C., and Smith, Gerald L.: Statistical Filter Theory 

For Time-Varying Systems. I.R.E. National Convention Record, 

Part 4, 1959j or I.R.E. Transactions on Automatic Control, vol. 

AC- 4, no. 2, Nov. 1959 , pp. 7^79- 

7 . Shinbrot, Marvin, and Carpenter, Grace C.: An Analysis of the 

Optimization of a Beam Rider Missile System. NACA TN 4l45, 1958* 

8 . Mersman, W. A.: Evaluation of an Integral Occurring in Servomechanism 

Theory. Pacific Jour. Math., vol. 2, 1952, pp. 627 - 632 . 












- 






’ 








- 











































ro > 


A. Amplitude random 
Interval l random 

P(t) = (l/l) exp(-t/t) 



V'c (w) = 


X(X S ) 2 


7T ( a; 2 + X 2 ) 

X = Average switching rate 


B. Same as A except 

..2 2 
x s = x s = 0 T 


XsU) 


4 ) *“ 


a n 


X a. 
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7 r (cu 2 + X 2 ) 


X = Average switching rate 
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at end of each 
interval. 
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A. 
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Figure 1.- Example signal inputs. 
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(a) Error components. 

Figure 2.- Optimum performance of time -invariant systems and time-varying 
homing systems with stationary inputs. 
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(b) Total error. 


Figure 2.- Continued 
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(c) Acceleration components 
Figure 2.- Continued. 
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(d) Total acceleration. 


Figure 2 


Continued 
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(e) Total error as a function of total acceleration. 
Figure 2.- Concluded. 
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Figure 3 • - Optimum impulse responses 
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Figure 
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( a) Error . 


4.- Optimum performance for nonstationary step acceleration signal. 
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(t) Acceleration. 
Figure Continued. 
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(c) Total error as a function of total acceleration. 


Figure 4 .- Concluded 
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( a) Error component s 


Ei and E 2 . 


Figure 5 


Off-design performance of system optimized for stationary 
signal, actual signal nonstationary. 
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Figure 5«~ Continued. 
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(c) Acceleration components A x and A s . 
Figure 5 * _ 


Concluded 
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Figure 6 



V 

(a) Error component Ei. 

Off -design performance of system optimized for nonstationary 
signal, acutal signal stationary. 
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(t>) Acceleration component 
Figure 6.- Concluded. 
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Figure 7.- Effect of type of signal input for example case. 
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